   
	 tsset id year	 
	 xtbalance , range(1995 2018)
	 
	 	 
	 local x "LnCO2 LnEF LnTRA LnPGDP LnIS FI LnGI"
	 tabstat `x', s(mean p50 max min sd skewness kurtosis n) format(%6.4f) c(s)
	 
   	
		set seed 13379 
		 xtthres LnCO2 LnIS FI LnGI, thres(LnTRA) dthres(LnPGDP)  min(10) bs1(500) bs2(500) bs3(300)
		set seed 13791
		 xtthres LnEF LnIS FI LnGI, thres(LnTRA) dthres(LnPGDP)  min(10) bs1(500) bs2(500) bs3(300)
		 
		 		 
		xttr_graph  
       xttr_graph, m(22)  
       xttr_graph, m(21)
		 
//
	   dis e(rhat22)
	   dis e(rhat1)
       global q1 = min(e(rhat22),e(rhat1))   
	   global q2 = max(e(rhat22),e(rhat1))
	   dis "$q1" 
	   dis "$q2"  
	   dropvars d7 d8 LnPGDP_x_*  
	   gen d7 = (LnTRA<=$q1)               
	   gen d8 = (LnTRA> $q2)
	   gen LnPGDP_x_LnTRA1 = LnPGDP*d7             
	   gen LnPGDP_x_LnTRA3 = LnPGDP*d8
	   local x "LnIS FI LnGI"         
	   xtreg LnCO2 `x' LnPGDP_x_LnTRA* pgdp, fe        
	   est store fe
	   xtreg LnCO2 `x' LnPGDP_x_LnTRA* pgdp, fe robust 
	   est store fe_robust
	   local m "fe fe_robust"
	   esttab `m', mtitle(`m') nogap s(r2 r2_w N F) 
	          star(* 0.1 ** 0.05 *** 0.01)          
			  order(`x' LnPGDP_x_LnTRA1 LnPGDP LnPGDP_x_LnTRA3 )
			  
			  
//
	   dis e(rhat22)
	   dis e(rhat1)
       global q1 = min(e(rhat22),e(rhat1))   
	   global q2 = max(e(rhat22),e(rhat1))
	   dis "$q1" 
	   dis "$q2"  
	   dropvars d7 d8 LnPGDP_x_*  
	   gen d7 = (LnTRA<=$q1)               
	   gen d8 = (LnTRA> $q2)
	   gen LnPGDP_x_LnTRA1 = LnPGDP*d7             
	   gen LnPGDP_x_LnTRA3 = LnPGDP*d8
	   local x "LnIS FI LnGI"         
	   xtreg LnEF `x' LnPGDP_x_LnTRA* pgdp, fe        
	   est store fe
	   xtreg LnEF `x' LnPGDP_x_LnTRA* pgdp, fe robust 
	   est store fe_robust
	   local m "fe fe_robust"
	   esttab `m', mtitle(`m') nogap s(r2 r2_w N F) 
	          star(* 0.1 ** 0.05 *** 0.01)          
			  order(`x' LnPGDP_x_LnTRA1 LnPGDP LnPGDP_x_LnTRA3 )
			  
 //
		 dis e(rhat1)
       global q1 = e(rhat1)  
	   dis "$q1"  
	   dropvars d5 d6 LnPGDP_x_*  
	   gen d5 = (LnTRA<=$q1)                 
	   gen d6 = (LnTRA> $q1)
	   gen LnPGDP_x_LnTRA1 = LnPGDP*d5              
	   gen LnPGDP_x_LnTRA3 = LnPGDP*d6
	   local x "LnIS FI LnGI"         
	   xtreg LnEF `x' LnPGDP_x_LnTRA* LnPGDP, fe       
	   est store fe
	   xtreg LnEF `x' LnPGDP_x_LnTRA* LnPGDP, fe robust 
	   est store fe_robust
	   local m "fe fe_robust"
	   esttab `m', mtitle(`m') nogap s(r2 r2_w N F)
	          star(* 0.1 ** 0.05 *** 0.01)          
			  order(`x' LnPGDP_x_LnTRA1 LnPGDP LnPGDP_x_LnTRA3 )
			  
			  
 //
		 dis e(rhat1)
       global q1 = e(rhat1)  
	   dis "$q1"  
	   dropvars d5 d6 LnPGDP_x_*  
	   gen d5 = (LnTRA<=$q1)                 
	   gen d6 = (LnTRA> $q1)
	   gen LnPGDP_x_LnTRA1 = LnPGDP*d5              
	   gen LnPGDP_x_LnTRA3 = LnPGDP*d6
	   local x "LnIS FI LnGI"         
	   xtreg LnCO2 `x' LnPGDP_x_LnTRA* LnPGDP, fe       
	   est store fe
	   xtreg LnCO2 `x' LnPGDP_x_LnTRA* LnPGDP, fe robust 
	   est store fe_robust
	   local m "fe fe_robust"
	   esttab `m', mtitle(`m') nogap s(r2 r2_w N F)
	          star(* 0.1 ** 0.05 *** 0.01)          
			  order(`x' LnPGDP_x_LnTRA1 LnPGDP LnPGDP_x_LnTRA3 )	  
			  

